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Abstract.  Far-infrared  spectroscopy  of  a  small  sample  of  IR-bright 
galaxies  taken  with  the  Infrared  Space  Observatory  Long  Wavelength 
Spectrometer  has  revealed  a  dramatic  progression  extending  from  strong 
fine-structure  line  emission  from  photoionized  and  photodissociated  gas 
in  the  starburst  galaxy  Arp  299  to  faint  [C  II]  158  pm  line  emission  and  ab¬ 
sorption  in  lines  of  OH,  H2O,  CH,  and  [O  I]  in  the  ultraluminous  galaxy 
Arp  220.  The  progression  to  weaker  emission  lines  is  accompanied  by 
lower  excitation  and  does  not  appear  to  be  due  to  far-infrared  extinction 
or  density  effects.  Although  aging  of  short  duration  starbursts  with  a  rel¬ 
atively  small  range  in  age  might  explain  much  of  the  sequence,  the  spectra 
of  Arp  220  and  other  ultraluminous  galaxies  may  be  more  reasonably  ex¬ 
plained  by  high  ionization  parameter,  dust-bounded  nebulae.  To  test  this 
hypothesis,  we  ran  photoionization  models  of  mid-  and  far-infrared  lines 
for  power  law  (index  -1.5)  and  starburst  cases  with  moderate  (500  cm“^) 
and  high  (10®  cm  densities.  Indeed,  we  find  that  for  dust-bounded 
nebulae  and  photoionization  by  a  power  law  source,  the  relative  emission 
line  strengths  and  the  high-to-low  ionization  infrared  line  ratios  decrease 
with  ionization  parameter.  For  starburst  ionization  sources  and  dust- 
bounded  nebulae,  the  [N  HI]57/L/m/[N  lV\122pm  line  ratio  is  found  to 
drop  with  ionization  parameter. 
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1.  Background 


Prior  to  the  launch  of  the  European  Space  Agency’s  Infrared  Space  Observatory 
(ISO)  in  1995,  Voit  (1992)  showed  how  mid-  and  far-infrared  fine-structure  lines 
could  be  used  to  constrain  the  electron  densities,  extinction,  and  shape  and 
ionization  parameters  of  the  central  ionizing  sources  in  ultraluminous  infrared 
galaxies  (ULIRGs).  Moreover,  the  ground-based  work  of  Roche  et  al.  (1991) 
showed  that  the  mid-IR  spectra  of  the  nuclei  of  galaxies  could  be  placed  into 
three  classes:  those  with  aromatic  feature  emission,  featureless,  and  those  with 
silicate  absorption,  typically  associated  with  optically  identified  H  II  region, 
Seyfert  1,  and  Seyfert  2  nuclei  respectively.  Building  on  this  early  work  on 
optically  selected  starburst  and  AGN  galaxies,  Genzel  et  al.  (1998)  constructed 
a  diagnostic  diagram  of  the  ratio  of  high-to-low  ionization  fine-structure  lines  vs. 
the  equivalent  width  of  the  7.7  nm  aromatic  feature  emission  based  on  which  they 
concluded  that  70  -  80%  of  ULIRGs  are  powered  predominantly  by  starbursts 
and  20  -  30%  are  powered  by  a  central  AGN.  They  attributed  the  weakness  of 
the  mid-infrared  fine-structure  lines  relative  to  the  infrared  luminosity  to  the 
effects  of  extinction. 

Far-infrared  spectroscopy  of  a  small  sample  of  IR-bright  and  ultraluminous 
galaxies  taken  with  the  ISO  Long  Wavelength  Spectrometer  (LWS)  has  revealed 
a  dramatic  progression  extending  from  strong  fine-structure  line  emission  from 
photoionized  and  photodissociated  gas  in  the  starburst  galaxy  Arp  299  (Satya- 
pal  et  al.  2001)  to  faint  [G  II]  158  nm  line  emission  and  absorption  in  lines  of  OH, 
H2O,  GH,  and  [O  I]  in  the  ULIRG  Arp  220  (Fischer  et  al.  1999).  With  the  pro¬ 
gression  towards  weak  emission  line  strengths,  the  sequence  shows  a  trend  toward 
lines  with  lower  excitation  potentials  in  the  ratios  [N  III]57 /:/m/[N  II]122  nm  and 
[O  III] 52  /:/m/[N  III] 57  fam.  No  FIR  fine-structure  line  emission  from  species  with 
excitation  potentials  greater  than  13.6  eV  were  detected  in  Arp  220.  These 
trends  do  not  appear  to  be  primarily  due  to  either  differential  far-infrared 
extinction  between  52  and  88/:/m  or  density  effects,  since  the  temperature- 
insensitive  [O  III]52/:/m/[0  III]88/:/m  line  ratio  does  not  show  a  trend  with  the 
ratio  [O  III]88/:/m/FIR  ratio  as  it  would  in  either  of  these  cases  and  all  of  the 
measured  [O  III]  line  ratios  fall  within  the  range  0.6  -  1.2,  consistent  with  elec¬ 
tron  densities  between  100  -  500  cm“^  (see  Fischer  et  al.  for  details).  Voit  (1992) 
discussed  the  possibility  that  even  the  mid-  and  far-infrared  fine-structure  lines 
would  be  weak  in  ULIRGs  if  they  are  formed  in  high  ionization  parameter  re¬ 
gions.  In  such  regions  with  high  ratios  of  ionizing  photon  to  electron  densities, 
UV  photons  are  preferentially  absorbed  by  dust  rather  than  gas,  due  to  the  high 
column  densities  of  ionized  gas  in  such  regions.  Bottorff  et  al.  (1998)  found  that 
the  L/^/H^  ratios  in  such  dust-bounded  nebulae  are  greater  than  100  for  ioniza¬ 
tion  parameters  greater  than  10“^,  density  =  100  cm“^,  and  stellar  temperatures 
between  30,000  -  50,000  K. 

Here  we  present  photoionization  models  with  starburst  and  power  law  ion¬ 
ization  sources  to  predict  the  strengths  of  the  fine-structure  lines  in  dust-bounded 
nebulae  and  to  compare  them  with  the  LWS  spectra  of  the  infrared-bright  galax¬ 
ies.  Due  to  the  weakness  of  the  infrared  fine-structure  lines  from  photoionized 
gas  in  ULIRGs,  ISO  mid-  and  far-infrared  spectroscopy  produced  mostly  upper 
limits.  Gomparison  of  photoionization  models  of  dust-bounded  nebulae  with 
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Figure  1.  Far-infrared  (left)  and  mid-infrared  (right)  line  ratios  vs. 
line-to-luminosity  ratios  for  an  AGN-like  power  law  ionizing  source  for 
moderate  (500  cm“^)  and  high  (10®  cm“^)  densities  are  plotted  for 
ionization  parameters  U,  for  log  U  =  -3,  -2,  -1,  0,  and  0.5  (as  shown 
by  the  greyscale  bar,  right).  See  text  for  the  details  of  the  modeling. 


spectra  from  future  space  missions  such  as  SIRTF  and  Herschel  will  help  to 
determine  the  conditions  in  the  photoionized  regions  of  these  galaxies. 

2.  Photoionization  modeling 

The  photoionization  modeling  was  done  using  CLOUDY  94.01  (Ferland  et  al. 
1998)  for  central  power  law  (Figure  1)  and  instantaneous  starburst  (Figure  2) 
ionization  sources.  For  power  law  models,  the  “table”  power  law  option  in 
CLOUDY  was  used.  This  option  produces  a  continuum  with  oc  that 
is  well  behaved  at  both  high  and  low  energy  limits  (10“®  -  10®  Rydbergs  in 
CLOUDY).  An  index  a  =  -1.5  was  used  for  the  mid-range  (10/L/m  -  50  keV), 
while  the  default  indices  of  -1-2.5  and  -2  were  used  for  the  low  and  high  ranges, 
respectively.  For  the  starburst  models,  we  used  the  instantaneous,  Salpeter  IMF, 
3  and  5  Myr  aged  burst  models  of  Leitherer  et  al.  (1999)  with  solar  metallicity 
and  standard  mass  loss.  H  II  region  abundances  were  used  in  CLOUDY. 

With  a  central  ionizing  source,  the  ionization  parameter  U,  defined  as  the 
ratio  of  ionizing  photons  to  hydrogen  atoms  at  the  inner  face  of  the  cloud,  is 
equal  to  Q/47rr^nc,  where  Q  is  the  central  Lyman  continuum  rate,  n  is  the 
density,  and  r  is  the  inner  radius  of  the  cloud.  It  was  varied  from  10“®  —  10^  by 
setting  Q  =  4.5  x  10®^  sec“^  and  varying  the  inner  radius  from  30  -  1600  pc. 

The  power  law  models  predict  a  decrease  in  the  [N  III]57/L/m/[N  II]122/L/m 
and  [Ne  V]14/L/m/[Ne  III]15/L/m  line  ratios  in  the  high  U  cases.  For  the  starburst 
models,  only  the  far-infrared  line  ratio  decreased  at  high  U,  over  the  range  of 
parameters  explored.  In  both  cases  the  line-to-luminosity  ratios  drop  at  high 
U,  as  expected.  These  effects  make  it  difficult  to  distinguish  between  AGN  and 
starbursts.  As  Figures  1  &  2  show,  high  densities  also  produce  both  lower  values 
of  the  line  ratios  and  lower  line-to-luminosity  ratios.  Older  starbursts  can  also 
produce  these  effects,  but  are  unlikely  to  power  ULIRGs  (Satyapal  et  al.  2001). 
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Figure  2.  Same  as  Figure  1  (except  for  log  scale  in  the  line  ratios), 
for  ionization  by  instantaneous  starbursts  of  ages  3  and  5  Myrs. 


3.  The  effects  of  high  ionization  parameters  on  photodissociation 
region  diagnostics 

Some  of  the  other  infrared  characteristics  of  ULIRGs  are  warm  60/100 /L/m  col¬ 
ors,  low  aromatic  feature-to-luminosity  and  low  [C  II]158/L/m  line-to-luminosity 
ratios,  although  their  aromatic  feature-to- [C  II]  ratios  are  normal  (Dudley  et  al. 
2001;  Luhman  et  al.  2001).  These  characteristics  may  be  the  result  of  the  ef¬ 
fects  of  grains  in  ionized  regions.  For  example,  Boulanger  et  al.  (1988)  find  that 
60/100 /L/m  colors  are  higher  and  small  grain  emission  lower,  as  traced  by  the 
IRAS  12  nm  flux,  within  the  California  nebula  H  II  region  than  outside  of  it. 
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